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Abstract 

The processes Y>(25) — > -jXcj, XcJ pp (J = 0, 1, 2) are studied using a sample of 14 x 10 6 Y>(25) 
decays collected with the Beijing Spectrometer at the Beijing Electron-Positron Collider. Very clear 
XcO, Xci and \c2 signals are observed, and the branching fractions B{\ c j — ► PP) (J = 0, 1,2) are 
determined to be (27.1±|;|±4.7) x 10~ 5 , (5.7^+0.9) x 10" 5 , and (6.5l| 4 ±1.0) x 10~ 5 , respectively, 
where the first errors are statistical and the second are systematic. 

PACS numbers: 13.25.Gv, 14.40. Gx, 12.38. Qk 
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I. INTRODUCTION 



Hadronic decay rates of P-wave quarkonium states provide good tests of quantum chromo- 
dynamics (QCD). The decays XcJ — ► PP have been calculated using different models [1, 2], 
and recently, the decay branching fractions of XcJ —> baryon and anti-baryon pairs were 
calculated including the contribution of the color-octet fock (COM) states [3]. Using the 
Xcj — ► PP branching fractions as input to determine the matrix element, the partial widths 
°f Xcj —> AA are predicted to be about half of those of Xcj —> PP, for J = 1 and 2. However, 
recent measurements of Xcj — * AA [4] together with the branching fractions of Xcj PP 
from pp annihilation experiments [5-7] and from a measurement from ip{2S) decays [8] seem 
to contradict this prediction. An improved measurement of the Xcj — ► PP branching fraction 
with the same data sample that was used for Xcj AA measurement will yield a more 
consistent measurement of this fraction. 

The measurements of B(xcj — > pp) have been performed in e + e~ collision experiments, 
where the Xcj are produced in ip(2S) radiative decays, and in pp annihilation experiments, 
where Xcj are formed directly. Although the precision in these experiments is limited, 
results from pp annihilation experiments seem systematically higher than those obtained in 
e + e~ — > ip{2S) experiments, as shown in Table I. This led to a global fit based on results 
from both e + e~ and pp annihilations [9]. 



TABLE I: B{xcj — ► pp) results obtained by different experiments. 





Channel 


Experimental technique 






ip(2S) -> 7x c j -> 1PP 


PP -> XcJ -> jJ/ip 






15.9±4.3±5.3 [8] 


48±i±?l [5] 
41 ± 3tl 6 [6] 


B(xci 


— ► pp)(xW~ 5 ) 


4.2±2.2±2.8 [8] 


7.8±1.0±1.1 [7] 


B{ X c2 




5.8±3.1±3.2 [8] 


9.1±0.8±1.4 [7] 



The above e + e~ annihilation results [8] were obtained from a sample of 3.79 x 10 6 ip(2S) 
events collected with the Beijing Spectrometer (BESI) detector [10] at the Beijing Electron- 
positron Collider (BEPC) storage ring running at the energy of the ip(2S). Here we report on 
a result obtained with a sample of (14 ±0.6) x 10 6 ip(2S) events collected with the upgraded 
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BESII detector [11]. In BESII, a 12-layer vertex chamber (VTC) surrounding the beam 
pipe provides trigger information. A forty-layer main drift chamber (MDC), located radially 
outside the VTC, provides trajectory and energy loss (dE/dx) information for charged tracks 
over 85% of the total solid angle. The momentum resolution is a p /p = 0.018a/1 + p 2 (p in 
GeV/c), and the dE/dx resolution for hadron tracks is ~ 8%. An array of 48 scintillation 
counters surrounding the MDC measures the time-of-flight (TOF) of charged tracks with 
a resolution of ~ 200 ps for hadrons. Radially outside the TOF system is a 12 radiation 
length, lead-gas barrel shower counter (BSC). This measures the energies of electrons and 
photons over ~ 80% of the total solid angle with an energy resolution of a E /E = 21%/y/E 
(E in GeV). Outside of the solenoidal coil, which provides a 0.4 Tesla magnetic field over 
the tracking volume, is an iron flux return that is instrumented with three double layers of 
counters that identify muons of momentum greater than 0.5 GeV/c. 

A Monte Carlo (MC) simulation is used for the determination of mass resolution and 
detection efficiency. The angular distribution of ip{2S) — > 7X is simulated assuming a pure 
El transition, namely 1 + cos 2 9, 1 — | cos 2 9, and 1 + ^ cos 2 9 for the spin-parity of the 
resonance X being J p = + , 1 + and 2 + , respectively, and X decaying to pp is simulated 
according to phase space. A Geant3 based package, SIMBES, is used for the simulation of 
detector response, where the interactions of the produced particles with the detector material 
are simulated and detailed consideration of the detector performance (such as dead electronic 
channels) is included. Reasonable agreement between data and Monte Carlo simulation has 
been observed in various channels tested, including Bhabha, e + e~ — > fi + fi~, J/ift — > pp, and 
ip(2S) -> 7T+7T- J/tjj, J/ijj -> g+f (£ = eoTfi). 

II. EVENT SELECTION 

To select 4>(2S) — > 'JXcj, Xcj —> pp (J = 0,1,2) candidates, events with at least one 
photon and two charged tracks are required. A neutral cluster in the BSC is considered to 
be a photon candidate when the angle between the nearest charged track and the cluster in 
the xy plane is greater than 15°, the first cell hit is in the beginning 6 radiation lengths, and 
the angle between the cluster development direction in the BSC and the photon emission 
direction in xy plane is less than 37°. 

A likelihood method is used for discriminating pion, kaon, proton, and antiproton tracks. 
For each charged track, an estimator is defined as W l = ^751, P l = Yij Pj{ x j)i where P % 

i 
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is the probability under the hypothesis of being type i, i — 7r, K, and p or p, and Pj(xj) 
is the probability density for the hypothesis of type i, associated with the discriminating 
variable Xj. Discriminating variables used for each charged track are the time of flight in 
the TOF and the energy loss of the track in the MDC. By definition, pion, kaon, proton and 
antiproton tracks have corresponding W l value near one. 

For the decay channel of interest, the candidate events are required to satisfy the following 
selection criteria: 

1. There are two oppositely charged tracks in the MDC with each track having a good 
helix fit and | cos#| < 0.75, where 6 is the polar angle of the track; 

2. At least one charged track is identified either as a proton or an anti-proton with 
W p > 0.7 or W p > 0.7; 

3. There is at least one photon candidate. In the case of multiple photon candidates, the 
one with the largest BSC energy is chosen as the photon radiated from the ip(2S); 

4. The x 2 probability of the four-constraint kinematic fit is required to be greater than 
1%. 

Figures la and lb show distributions of W p and W p after all other requirements have 
been applied. It can be seen that the W p ^ > 0.7 requirement rejects most of the tc and K 
background, while retaining high efficiency 

In order to reduce backgrounds from e + e~ or tracks in J/ip decays (ip(2S) — > 

XJ/ip), the BSC energy of the positively charged track (ESC P ) is required to be less than 0.7 
GeV, and the two charged tracks must satisfy the muon veto requirement, that N^ + N^ < 
6, where N htt is the number of mu-counter layers with matched hits and ranges from to 3, 
indicating not a muon (0), a weak (1), medium (2), or strongly (3) identified muon track [12]. 
Distributions of ESC P and N^ lt + N^ lt are shown in Figures lc and Id. 

After the above selection, the invariant mass of the proton and anti-proton is shown in 
Fig. 2a, where clear Xco, Xa and Xc2 signals can be seen. The large peak near the mass of 
the ip(2S) is due to ip(2S) — > pp with a fake photon reconstructed. 

The same analysis is performed on a MC sample with 14 M inclusive ip(2S) decays 
generated with Lundcharm [13]. It is found that the remaining backgrounds are mainly 
from ip{2S) — > n°pp with tt° — > 27, and they contribute a smooth part to the pp invariant 
mass distribution. 
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FIG. 1: (a) W p ,(b) W p ,(c) ESC P , and (d) N^ + Ng* distributions after all the other requirements 
have been applied, as described in the text. The histograms are for Monte Carlo simulated ip(2S) — ► 
7XcO> XcO PP events and dots with error bars are for data with pp invariant mass within the Xcj 
signal region. For (a) and (b), data and MC simulation are normalized to the last bin, and for (c) 
and (d), data and MC simulation are normalized to ESC P < 0.7 and N f p nt + Np lt < 6, respectively. 

III. FIT TO THE pp INVARIANT MASS SPECTRUM 

The pp mass spectrum of the final selected events with pp mass between 3.26 and 
3.64 GeV/c 2 is fitted with three Breit-Wi gner resonances smeared by Gaussian mass reso- 
lution functions together with a second order polynomial background using the unbinned 
maximum likelihood method. In the fit, the mass resolutions are fixed to values from Monte 
Carlo simulation (6.84, 6.59 and 6.17 MeV/c 2 for Xco, Xci an d Xc2, respectively), and the 
widths of Xd and Xc2 are fixed to 0.92 and 2.08 MeV/c 2 coming from the PDG2002 [14]. 
The fit, shown in Fig. 2a, yields a likelihood probability of 79%, total numbers of events of 
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FIG. 2: Breit-Wigner fit to pp mass distribution for selected ip{2S) — > 7pp events (a) in data and 
(b) in Monte Carlo simulation, with relative branching fractions fixed to data. 

89.5l{3, 18.2I4J, and I4.3I4; 2 ., and statistical significances of 10.6<j, 4.6a, and 3.7a for Xco, 
Xd and Xc2, respectively. The fitted masses are 3414.3 ± 1.6, 3513.0 ± 2.1, and 3549. ll^Q 
MeV/c 2 , respectively, and agree well with the world averages [14], and the width of the 
Xco is determined to be I2.8I4J MeV/c 2 , in good agreement with results from other experi- 
ments [14]. 

The same fit is applied to the Monte Carlo sample, which is about 180 times larger than 
the data sample and uses the same relative branching fractions between the three Xcj states 
as determined from the data. The background fraction is estimated using the Xco mass 
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region for data. The fit yields the efficiencies for each channel as e Xc0 = (27.49 ± 0.30)%, 
e Xcl = (27.42 ± 0.56)%, and e Xc2 = (23.26 ± 0.50)%, where the errors are due to the limited 
statistics of the Monte Carlo samples. 

IV. SYSTEMATIC ERRORS 

Systematic errors of the measured branching fractions come from the efficiencies of photon 
identification, particle identification, the kinematic fit, etc. 

A. Photon identification 

To investigate the systematic error associated with the fake photon misidentification, the 
fake photon multiplicity distributions and the energy spectra in both data and Monte Carlo 
sample are checked with ip(2S) — >■ pp events. In this channel, it is found that the Monte 
Carlo simulates a little more fake photons than data. 

The selection of the photon candidate with the largest BSC energy has an efficiency of 
(92.55 ± 0.52)% for data and (94.49 ± 0.22)% for Monte Carlo in ip(2S) -> jXc2 -> 1PP, as 
determined from the comparison of energy distributions between the real photons and the 
fake ones; the difference is (2.1 ± 0.6)%. For Xci and Xco, the photon is more energetic, so 
the efficiency of selecting the largest BSC energy cluster is higher. The systematic error on 
the photon identification is taken as 2.7% for Xc2, as well as for Xco and xa- 

B. Photon detection efficiency 

The detection efficiency of low energy photons is studied with J/ip — > 7r + 7r~7r° events by 
requiring only one photon in the kinematic fit and examining the detector response in the 
direction of emission of the second photon. The efficiency from the Monte Carlo simulation 
agrees with data within 2% in the full energy range. This is taken as the systematic error 
of the photon detection efficiency. 

C. Particle identification 

The systematic error from particle identification is studied using p (or p) samples from 
ip(2S) — > ty + ty^J/iIj, J ftp — > pp and ip(2S) — > pp. Since only one track is required to be 
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identified, the efficiency is very high. The correction factors for efficiencies for Xco, Xa and 
Xc2 are found to be 1.020±0.006, 1.025±0.006 and 1.028±0.006 respectively. 

D. Kinematic fit 

The systematic error associated with the kinematic fit is caused by differences between 
the measurements of the momenta and the error matrices of the track fitting of the charged 
tracks and the measurement of the energy and direction of the neutral track for the data 
and the simulation sample. The effect is studied for charged tracks and neutral tracks 
separately. The systematic error due to the charged tracks was checked using ip(2S) — > pp 
events, which can be selected easily without using any kinematic fit. By comparing the 
numbers of the events before and after the kinematic fit, the efficiencies for x 2 probability 
> 1% are measured to be (85.34 ± 1.64)% for data and (88.17 ± 0.56)% for Monte Carlo 
simulation, respectively. This results in a correction factor of (0.968 ± 0.020)% for the 
efficiency of this specific channel. 

The uncertainty due to the measurement of the neutral track parameters is taken from 
Ref. [4] based on a study of i/)(2S) — > 7Xcj, Xcj — > Ti + ^pp- A systematic error of 4.2% is 
quoted for all the channels. 

E. Muon veto and BSC energy requirements 

The efficiencies of the muon veto and BSC energy requirements are studied with ip(2S) — > 
pp events. In order to get a clean sample of pp without using these requirements, a four- 
constraint kinematic fit to pp is done, and the x 2 probability of the fit is required to be greater 
than 1%. The efficiencies of the muon veto for both data and Monte Carlo simulation 
are found to be 100%, while those of the BSC energy requirement are (98.94 ± 0.40)% 
and (99.14 ± 0.18)% for data and Monte Carlo, respectively, resulting in a difference of 
(0.20 ± 0.44)%. This difference together with the error is taken as the systematic error of 
the BSC energy requirement. 
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F. Angular distribution 



For the radiative decay 4>(2S) — > 'JXcj, the general form of the angular distribution is 
W(cos9) = 1 + A cos 2 9, where 9 is the angle between the beam direction and the outgoing 
photon and \A\ < 1 [15]. The value of A can be unambiguously predicted only for spin 
J = 0, where A — 1. By comparing the multipole coefficients measured by Crystal Ball [16] 
with those used under the assumption of pure El transition, differences of 3.5% and 5.5% 
are found, which will be taken as systematic errors of the angular distributions for J = 1 
and 2 respectively. 

G. Breit-Wigner fit 

To determine the fit systematic errors, different background shapes and fit ranges are 
used in the fit of the pp invariant mass distribution. After changing the background shape 
from a second order to a first order polynomial and varying the fitting range around the one 
used in the fit, the uncertainties due to fit are determined to be 11.5%, 8.4% and 7.3% for 
X c o, Xd, and Xc2, respectively. 

H. Monte Carlo determined mass resolution 

The systematic error due to the use of the MC determined mass resolution is studied with 
ip(2S) — > it + it~J/vI), J/ip — > pp. The non-Gaussian tails are found to be (5.3 ± 4.4)% and 
(5.3 ±1.9)% for data and Monte Carlo respectively, which indicates good agreement between 
data and Monte Carlo simulation. The uncertainty of the comparison, 4.8%, is taken as the 
systematic error due to the MC determined mass resolution. 

I. Other systematic errors 

The results reported here are based on a data sample corresponding to a total number of 
i/)(2S) decays, N^s), of (14.0±0.6) x 10 6 , as determined from inclusive hadronic events [17]. 
The uncertainty of the number of i/)(2S) events, 4%, is determined from the uncertainty in 
selecting the inclusive hadrons. 

The difference of MDC tracking efficiencies between data and Monte Carlo for p and p 
may cause a systematic error of 1-2% for each track. Here 3% is taken as the systematic error 



10 



on the overall tracking efficiency. The trigger efficiency is around 100% with an uncertainty 
of 0.5%, as estimated from Bhabha and e + e~ — > events. The systematic errors on the 

branching fractions used are obtained from PDG [14] directly. 

J. Total systematic error 

Table II lists the systematic errors from all sources. Adding all errors in quadrature, the 
total errors are 17.3%, 15.4% and 15.6% for Xco, Xci and x c2 , respectively. The corresponding 
correction factors (/) for the Monte Carlo simulated efficiencies are 0.987, 0.992, and 0.995. 

TABLE II: Summary of systematic errors in percent. Numbers common for all channels are only 
listed once. 



Source 


XcO 


Xcl 


Xc2 


MC statistics 


1.1 


2.0 


2.2 


Photon ID. 




2.7 




Photon eff. 




2 




4C-fit (neutral) 




<4.2 




4C-fit (charged) 




2.0 




Particle I.D. 




0.6 




BW fit 


11.5 


8.4 


7.3 


Mass resolution 




4.8 




BSC energy cut 




0.6 




Angular distr. 





3.5 


5.5 


Number of tp(2S) 




4 




MDC Tracking 




3 




Trigger Efficiency 




0.5 




BmS) -+ 7Xcj) 


9.2 


8.3 


8.8 


Total Systematic error 


17.3 15.4 15.6 



11 



V. RESULTS 



The branching fraction of Xcj PP is calculated using 

n° bs /{e-f) 



B(Xcj -> PP) = 



N H2S ) ■ B[iP(2S) - 7Xcj] 
Using numbers listed in Table. Ill, one obtains 



&(Xco ^PP) 
B(xa -> PP") 

#(Xc2 -^PP~) 



(27.lt£;jj ± 4.7) x 10~ 5 , 
(5.7+^ ±0.9) x 10" 5 , 
(6.5^ ±1.0) x 10- 5 , 



where the first errors are statistical and the second are systematic. The measured branching 
fractions agree with corresponding world averages within errors [14]. 

TABLE III: Numbers used in branching fraction calculation and the final results. 



quantity 


XcO 


Xcl 


Xc2 


n obs 

e(%) 


89.5^ 
27.49±0.30 


18.2±t;| 
27.42±0.56 


14.3±f 7 
23.26±0.50 


f 


0.987 


0.992 


0.995 


^(25)(10 6 ) 
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B[^{2S) ^ 1X cj\(%) 


(8.7±0.8)% 


(8.4±0.7)% 


(6.8±0.6)% 


B(Xcj - VP) (10- 5 ) 
Rb 


± 4.7 5.7t\i ± 0.9 G.5±H ± 1.0 
1.73 ±0.63 4.56 ±2.34 5.08 ± 3.08 



The relative branching fraction of Xcj AA to Xcj PP is found with the following 
formula: 

3 n obs If - 

PP I c pp 

These results are also shown in Table III, by using the numbers in Table III of this paper 
and those in Table II of Ref. [4], with the common errors in B(x c j AA) and B(x c j PP) 
canceled out. The measurements confirm the enhancement of Xcj — * AA relative to Xcj 
pp, as compared with the COM calculation [3]. 
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The actual measured quantities in this analysis is the branching fractions of i/)(2S) 
IXcj -> iPP, with 



Bty(2S) - 7XcJ - iPP] = B[^{2S) - 7Xcj] • £(x c j - pp) 

n obs /(e ■ f) 

using numbers in Table III, the results are 

B[^(2S) - 7Xco - 7PP] = (23.6+H ± 3-4) x 1(T 6 , 
B[i/>(2S) -> 7Xd - 7PPl = (4.8+l;£ ± 0.6) x lO" 6 , 
6^(25) -> 7Xc2 -> 7PP] = (4-4lii ± 0.6) x 10' 6 . 

VI. SUMMARY 

Decays Xcj - ► are observed using the BESII sample of 14 million ip(2S) events, and the 
corresponding branching fractions are determined. The measured values agree with previous 
experiments within errors [5-8]. The measurements confirm the enhancement of Xcj A A 
relative to Xcj — ► PP, as compared with the COM calculation [3] . 
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